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We develop a theoretical treatment of polarization spectroscopy and use it to make predictions about the 
general form of polarization spectra in the alkali atoms. Using our model, we generate theoretical spectra for 
the D2 transitions in Rb, Rb, and Cs. Experiments demonstrate that the model accurately reproduces 
spectra of transitions from the upper hyperfine level of the ground state only. Among these, the closed transition 
F — > F' — F + 1 dominates, with a steep gradient through line center ideally suited for use as a reference in 
laser locking. 
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I. INTRODUCTION 

Polarization spectroscopy is a sub-Doppler spectroscopic 
technique in which birefringence is induced in a medium by a 
circularly polarized pump beam, and interrogated with a coun- 
terpropagating weak probe beam It is closely related 

to saturation spectroscopy, and has proved a useful tool in ar- 
eas as diverse as combustion diagnostics |4], plasma charac- 
terization |5], and laser frequency stabilization |6, 7, 8]. 

In this work we develop a simple theoretical model of polar- 
ization spectroscopy based on a calculation of rate equations. 
Using this model, we generate spectra for the hyperfine tran- 
sitions of the D2 lines of rubidium and cesium, and show that 
for transitions from the upper hyperfine level of the ground 
state (n 2 S 1/2 (F =1+ 1/2) -> n 2 P 3/2 (F' = F+l, F, F—l, 
where 1 is the nuclear spin) our model gives good agreement 
with the observed spectra. For transitions from the lower hy- 
perfine level, the model breaks down, demonstrating that fac- 
tors not included in the model play a significant role. 

The layout of this paper is as follows. In section 2 we de- 
rive expressions needed to construct our model, present results 
for the time-dependence of the population of relevant hyper- 
fine states, and discuss how spectra are produced. Section 3 
describes the experimental apparatus used to test the model's 
predictions in rubidium and cesium, and section 4 contains 
experimental and theoretical spectra. Section 5 contains a dis- 
cussion of our results, and in section 6 we draw our conclu- 
sions. 



II. THEORY 

A. Optical Bloch and population change equations 

The density matrix is used to calculate the atom-light in- 
teraction. Consider a two level atom with ground state a, 
excited state b, and transition angular frequency u>o, inter- 
acting with laser light of angular frequency wl and detuning 
A = cjl — wq. The Optical Bloch equations are |9] 
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Here p aa and pbb are the ground and excited state proba- 
bilities, respectively, and the off-diagonal terms p a b and pba 
are coherences. (The tilde notation means that we are using 
"slow" variables - i.e. the evolution of the coherence at the 
laser frequency has been factored out separately). The Rabi 
frequency is denoted by Or and is proportional to the prod- 
uct of the electric field of the laser and the transition matrix 
element. The effect of spontaneous emission is to reduce the 
excited state population at a rate T = 1 /r, and the coherences 
at a rate T/2, where t is the excited state lifetime |9]. There 
are no collisional terms as the mean time between collisions in 
the room temperature vapor cells used in these experiments is 
far longer than the mean time spent in the laser beams. Similar 
equations are derived for many-level atoms for every Zeeman 
sub-state. 

It becomes much easier to solve these equations if the co- 
herences can be eliminated. To show that such an elimina- 
tion is valid for our system, we numerically solved the optical 
Bloch equations containing the coherences and the rate equa- 
tions for a few sample systems. For a short time the solutions 
differ; the rate equations are linear in time whereas the opti- 
cal Bloch equations give a quadratic time dependence for the 
evolution of populations. However, the solutions become in- 
distinguishable after a time of ~ 7r 11(11 . where t=26.24(4) 
ns for rubidium |11] and 30.57(7) ns for cesium Hal . As 
the typical time of flight of an atom through the laser beams 
(discussed below) is roughly two orders of magnitude longer 
(~ 5 x 10~ 6 s) it is an excellent approximation to solve the 
equations for the rate of change of populations. We there- 
fore assume that the coherences evolve sufficiently quickly 
that their steady state value can be used. 

With the coherences eliminated, we set the time derivatives 
to zero in the second and third equations of Q and substitute 
into the first equation, obtaining 
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This leaves rate equations for populations only. Recall- 
ing that the light intensity I is proportional to the square of 
the Rabi frequency, the saturation intensity I s &t is defined as 
I /hat = 2£!p/r 2 . The three terms on the right-hand side 
of (0 have the usual physical interpretation: the first repre- 
sents stimulated emission out of the excited state, which is 
proportional to the light intensity; the second absorption into 
the excited state, also proportional to the light intensity; and 
the third spontaneous emission out of the excited state, which 
is independent of the light intensity. The intensity-dependent 
terms have a Lorentzian line-shape, with full width at half 
maximum of T. Note that the usual rate equations of laser 
physics for populations with Einstein A and B coefficients as- 
sume broad-band radiation, and use the intensity per band- 
width (evaluated at line center); in contrast, equation (0 is 
valid for narrow -band radiation and is a function of the detun- 
ing. 

It is straightforward to generalize these equations for mul- 
tilevel systems. As there are three possible polarizations, and 
the AF = 0, ±1 selection rule applies, each excited state can 
decay into up to nine different ground states. The pump beam 
drives <r + transitions, and the Rabi frequencies and saturation 
intensities are calculated from the known line-strength coeffi- 
cients. 

The line strength for a transition between two states a and 
b is proportional to the square of the dipole matrix element, 
given 1 13] by 
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where L and S are respectively the orbital and spin angular 
momenta of the electron, J is their sum, F = X + J where X 
is the nuclear spin, and mp is the projection of F onto the z 
axis. The expressions in curly brackets are standard 6 J sym- 
bols encountered when re-coupling angular momenta, and the 
term with large round brackets is a 3 J symbol. These can 
be calculated using standard symbolic mathematical packages 
J3- The value of q denotes the change in z-component of the 
angular momentum in the transition between the two levels, 
with q = for tt transitions, and q = ±1 for a ± transitions. 

As an example, consider 87 Rb (X = 3/2), which has two 
hyperfine states in the ground term, F = 2, 1. If the pump 
laser is tuned in the vicinity of the resonances Si/ 2 (F = 

2) -> Pz/2(F' = 3, 2, 1) with the detuning, A, defined rel- 
ative to the closed transition Si/ 2 (F = 2) — » Pa/^F' — 

3) , the strongest transition will be 2 Si/ 2 (F — 2,mp = 
2) —> 2 P3/2(F' = 3, TOf' = 3). For this transition, the sat- 
uration intensity 7 sat has the value 1.6 mW/cm 2 . The reduced 



matrix element for the transition (55||r||5P) = 5.14ao, 
where ao is the Bohr radius, is calculated from the excited 
state lifetime. The line-strength ratios, which we define as 
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are always less than one as they are normalized relative to the 
strongest transition. 

To simplify our notation we write the population of a state 
as PF.m F , which is equal to the diagonal density-matrix ele- 
ment pF,mF.F,m F - The five F = 2, nip ground states obey 
the rate equations 
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where A' = A if F' = F + 1; A' = A + A 32 if F' = F; 
and A' = A + A 3i if F' = F - 1. Here hA 32 and hA 3 i 
are the excited state hyperfine intervals. The first row of Eq. 
contains the stimulated absorption and emission terms, the 
second row the spontaneous emission terms. 

The three F = l,mj? ground states obey the rate equations 
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(3) (6) 
Note that there are no stimulated terms as the large ground 
State splitting means that these transitions very far off reso- 
nance; these states can only increase their population by spon- 
taneous emission from the excited states. 

Finally, the excited state population rate equations are 
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A useful check before solving these equations is that the 
sum of all of the right-hand sides of these equations must 
be zero, because the total population in all states is constant. 
Note that there are 36, 24, or 48 Zeeman levels which have to 
be considered for the 85 Rb, 87 Rb, and 133 Cs systems, respec- 
tively; hence the simplification of the equations by eliminating 
the coherences is substantial. 
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B. Calculation of anisotropy 

The birefringence of the medium can be probed by a lin- 
early polarized probe beam, which can be decomposed into 
two beams of equal amplitude and opposite circular polariza- 
tion. For an isotropic medium, both circular components ex- 
perience the same refractive index and absorption. There is 
a difference in absorption coefficients, Aa, for an optically 
pumped medium: Aa = a + — a_, with a± being the ab- 
sorption coefficients of the circular components driving a 
transitions. Correspondingly, the incident linearly polarized 
beam exits the medium with a rotated elliptical polarization. 
Both the ellipticity and the angle of rotation are proportional 
to Aa |6]. In our experimental set-up (Section 3) we are sen- 
sitive to the rotation of polarization, providing the absorption 
coefficients are not too large. The signal has a characteristic 
dispersion spectrum (equation (9) in [6]). Therefore for each 
resonance we calculate one parameter, the line-center differ- 
ence in absorption. 

The experiments were performed with alkali metal atoms 
( 85 Rb, 87 Rb, 133 Cs) on the D2 transition. For an atom in the 
2 Si/ 2 ground state with nuclear spin I there are two values 
for the total angular momentum F, namely F = X ± 1/2. 
There are four values for the excited state angular momen- 
tum F', namely F' = 1 ± 3/2,1 ± 1/2. The hyperfine 
splitting of the ground states (3.0 GHz for 85 Rb, 6.8 GHz 
for 87 Rb, 9.2 GHz for 133 Cs) exceeds the room temperature 
Doppler width (~ 0.5 GHz), whereas the excited state hyper- 
fine splitting is less than the Doppler width. Therefore, the 
absorption spectrum consists of two isolated Doppler broad- 
ened absorption lines per isotope. Sub-Doppler absorption 
features are obtained by using counter-propagating pump and 
probe beams. There are three excited state resonances coupled 
to each ground state via electric dipole transitions, namely 
Si/ 2 (F) -> P3/2(F' = F+ 1,F,F - 1). In addition, one 
observes three cross-over resonances at frequencies halfway 
between each pair of conventional resonances fl5Tl . 

The calculation assumes that there is no excited state popu- 
lation initially, and that the ground state population is spread 
uniformly amongst the 2(2Z-fl) different \F, nip) levels. The 
pump beam has intensity / and drives er + transitions. This 
optical pumping process drives the population towards the 
largest possible value of nip, and the medium becomes op- 
tically anisotropic. Equations similar to Q are written down 
for each m p level, and the set of coupled equations are solved 
numerically. 

From these solutions, we obtain graphs of the population 
of each hyperfine state as a function of time. Figure^shows 
the time-evolution of these populations for a sample of 87 Rb 
atoms when A/27r=0, -267 MHz, and -424 MHz, correspond- 
ing to a laser on resonance with the F = 2^F' = 3,2,1 
transitions respectively. When the laser is resonant with the 
closed F = 2 — > F' — 3 transition, atoms are rapidly pumped 
into the extreme |2, 2) state, with over 50% of atoms in this 
state after <10 fj,s. The population of the other states is either 
constant or decreasing after ~ 1 /is. 

These figures clearly show the different optical pumping 
dynamics as a function of the laser detuning. A detuning 
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FIG. 1: Time-evolution of population. The population of the eight 
ground state sub-levels of 8r Rb are plotted as a function of time after 
the pump beam is turned on. The pump drives a + transitions and has 
an intensity of 0. 1 J sat • Vertical axes are identical to emphasize dif- 
ferences in population behavior, a) Laser tuned to F = 2 — > F' — 3 
transition. Within a few /j,s the initial isotropic population distribu- 
tion becomes highly anisotropic, as atoms are optically pumped into 
the |2, 2) state. As this transition is closed, no change occurs in the 
population of the F = 1 manifold, b) When the laser is tuned to the 
F — 2 — > F 1 — 2 transition the atoms are again optically pumped 
towards the |2, 2) state. However, as this is an open transition, a sig- 
nificant fraction of population accumulates in the F = 1 manifold. 
These states do not contribute to the medium's optical anisotropy. c) 
For the laser turned to the F = 2 — > F' = 1 transition even fewer 
atoms are pumped into the |2, 2) state. 
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A=0 (Figure^) is optimum for achieving the most optically 
anisotropic medium. For A/27r=-267 MHz (Figure [Q)) the 
medium becomes optically anisotropic, but to a lesser extent. 
This is due to the open nature of the transition, which allows 
atoms to accumulate in the F = 1 manifold, where they are 
too far from resonance to influence the medium's anisotropy. 
This effect is even more pronounced for A/27r=-424 MHz 
(Figure[l];), and the medium's optical anisotropy at this detun- 
ing is correspondingly reduced in comparison with the other 
detuning s. 

These population graphs allow us to predict that the polar- 
ization spectroscopy signal should be dominated by the closed 
transitions, S 1/2 (F = J + 1/2) -> P 3/2 {F' = F + 1) 
and S 1/2 {F =1-1/2) -► P 3/2 (F' = F - 1). For such 
transitions, selection rules forbid atoms in the excited state 
from falling into the other ground state; consequently, all of 
the ground state population ends up pumped into the extreme 
magnetic sublevels with nip = F (forX + 1/2, as discussed 
in the example above) ormj? = F and nip = F — 1 (for 
X — 1/2). The absorption coefficients ot± differ most for this 
extreme state; therefore, the anisotropy of the medium is max- 
imized. Note that the sign of the anisotropy generated by the 
1+1/2 transitions is expected to be opposite that of the lower 
hyperfine transitions, due to the lack of allowed er + transitions 
from the ground nip = F = X — 1/2 states. Figure [2] illus- 
trates the optical pumping process. 

After solving the coupled population rate equations it is 
possible to calculate the time-dependent anisotropy of the 
medium, A(t). This is defined as a sum over all of the ground 
itif states of the difference in the absorption coefficients for 
the components of the probe beam driving u ± transitions, tak- 
ing into account the ground and excited state populations: 

m F =+F 

-^(O = J! "^"if^'^f+i) (PF, mF - P'f>,w,f + i) 

171F— — F 



a) 




FIG. 2: Sign of anisotropy for closed transitions, a) For the closed 
F = 1+1/2^F' = F + 1 transition the population is optically 
pumped into the \F, thf = F) state, with a small fraction in the 
excited state. The line strength aF,m F ->F+i,m F + i is significantly 
larger than b) For the closed F = X — 1/2 — > 

F' = F — 1 transition the population is optically pumped into the 
\F, m,F = F) and \F, my = F — 1) states. There are no allowed 
<t + transitions for these states, whereas the line strength for the o~ 
transitions are finite. Consequently, the anisotropy of the medium 
has opposite sign relative to a). All levels are drawn as for an atom 
with 1 = 5/2; the structure of the other alkali atoms is similar. 
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For each isotope and each ground state F the anisotropy is 
calculated for the three frequencies corresponding to the res- 
onances F ->• (F' = F + 1, F, F - 1), i.e. A=0, -A 32 , and 
— A31. It is assumed that since the experiment uses counter 
propagating pump and probe beams only atoms with zero ve- 
locity along the axis of the beams contribute significantly. The 
details of the transverse motion of the atoms is outlined in the 
next sub-section. 



C. Transverse motion of atoms 

The experiment measures the average anisotropy of the 
medium. To simulate this, we average the time-dependent 
anisotropy, A(t), with a weighting function, 7i{t), which 
gives the distribution of times of flight transverse to the beam. 
For a circular beam of radius a the probability distribution, 



F(£), of having a path length £ in a uniform gas is 
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The probability distribution function, Q(t,£), for having a 
transit time t for a given £ is obtained from the Maxwell ve- 
locity distribution for a sample of atoms of mass m at temper- 
ature T and is 
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The distribution function H(t) is obtained thus: 



H(t)= / g(t,i)F(£)dt, 

J£=0 
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which can either be evaluated numerically or in closed form in 
terms of the complex error function. The average anisotropy 
A is calculated by averaging the time-dependent anisotropy 
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with the time of flight weighting function, i.e. 

A I A{t)H(t)dt. 



(12) 



This integral is calculated numerically at the three resonant 
frequencies of each D2 transition. 



D. Generating theoretical spectra 

To generate the predicted spectra, the three values of A 
from il 2i are multiplied by dispersion functions of the form 
x/(l + x 2 ), where x = 2A/T. As in the work of Yoshikawa 
et al. 1 16] the strength of cross-over features is assumed to 
be the average of the two associated resonances; hence cross- 
over features are calculated by multiplying this average by a 
dispersion function located halfway between each resonance. 

To account for spectral broadening produced by satura- 
tion effects (power broadening), we substituted a broadened 
linewidth 7 = kT MHz for the natural linewidth T in the 
dispersion function. We found that for 1 < k < 2, the ex- 
perimental data do not provide a tight constraint on the value 
of 7. Doppler broadening is incorporated by convolving the 
resulting spectra with a Gaussian of FWHM < 1 x 2w MHz, 
consistent with residual Doppler broadening due to the finite 
crossing angle between the probe and pump beams. Magni- 
tudes for both types of broadening were set by the experimen- 
tal conditions, as were the values for temperature (293K for 
rubidium, 273K for cesium), beam radii, and I/I aa ,t- 
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FIG. 3: Experimental layout. A thick glass slide picks off a fraction 
of the extended cavity diode laser (ECDL) light and splits it into 
two parallel beams. One beam acts as the probe, and the other as 
a (nearly) counterpropagating pump. The beams overlap inside a 
70 mm (50 mm) Rb (Cs) cell, which rests inside a long solenoid 
and below a coil that cancels the ambient laboratory field. In the 
cesium experiment two coils in the Helmholtz configuration replace 
the solenoid. A half-wave plate rotates the plane of polarization of 
the probe beam with respect to the polarizing beam splitter (PBS) 
cube axis; a quarter- wave plate makes the pump circularly polarized. 



TABLE I: Specifications of lasers used in this work. 



Atom 


A 


Horizontal r ( 1 /e 2 ) 


Vertical r (1/e 2 ) 


Rb 


780 nm 


0.65±0.01 mm 


0.59±0.01 mm 


Cs 


852 nm 


0.69±0.01 mm 


0.70±0.01 mm 



III. EXPERIMENTAL SETUP 



The layout of our experiment is similar to that described in 
J6J] and is shown in Figure[3] The probe beam's plane of polar- 
ization is set at ir/4 with respect to a polarizing beam splitter 
(PBS), which acts as an analyzer. The signal we record is 
the difference between the signals in each arm of the PBS; in 
the absence of a pump beam, the two arms will have equal 
intensities, and the difference will be zero. This technique 
produces polarization signals an order of magnitude larger 
than the conventional (single detector) method of polarization 
spectroscopy - an important advantage if the signal is being 
used for laser locking [6]. 

The experimental layouts for rubidium and cesium spec- 
troscopy are very similar, differing only in the length of the 
alkali vapor cell, equipment used to control the magnetic field 
along the axis of the cell, and waveplate type. Table 1 contains 
information on the spot radii (1/e 2 intensities) of the beams 
after the light has passed through a pair of anamorphic prisms. 

In both experiments the pump and probe beams are derived 
from the same extended cavity diode laser (ECDL). The diode 
for the rubidium (cesium) laser was a Sanyo DL-7 140-201 
(SDL-5401-G1). The crossing angle between probe and coun- 
terpropagating pump within the vapor cell is <3.0±0.2 mrad. 
Neutral density filters are used to vary pump and probe powers 
independently. A half-wave plate rotates the polarization of 
the probe relative to the axis of the PBS; a quarter-wave plate 



converts light in the pump beam to circular polarization. The 
two output beams from the PBS are focused onto photodiodes, 
which are connected to simple current-to-voltage circuits de- 
signed to output a voltage linearly proportional to the incident 
intensity. These voltages are then subtracted electronically to 
yield the polarization spectra. 

Since light of a given polarization may drive er + , er~ or 
it transitions depending on the external magnetic field, it is 
necessary to establish a "preferred" magnetic field direction 
along the vapor cell axis. In the rubidium experiment this is 
done by placing a room-temperature cell containing 85 Rb and 
8 'Rb inside a 300-turn solenoid of length 280 mm and diam- 
eter 26 mm. Numerical simulations showed that the magnetic 
field inside the solenoid is uniform to within 0.2% over the 
length of the 70 mm cell. 

In the cesium experiment the cell is partially submerged 
in an ice bath, and two 300 x 300 mm square coils in the 
Helmholtz configuration generate an axial magnetic field uni- 
form to 1%. The ice bath is needed because our model as- 
sumes an optically thin medium, but at 23 °C absorption ex- 
ceeds 90% for a 50 mm Cs vapor cell (compared to a maxi- 
mum of 30% for Rb in a 70 mm cell). In the ice bath this is 
reduced to 60% (50%) for transitions from the upper (lower) 
hyperfine level of the ground state. In both experiments we 
also cancel the (primarily vertical) ambient laboratory mag- 
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netic field with a 245 mm diameter coil mounted above the 
cell. Any inhomogeneity due to the finite diameter of this coil 
makes a negligible contribution to the total field when added 
in quadrature with the axial field. The magnitude of the axial 
field is set to a value just below the point where Zeeman split- 
ting of the hyperfine levels begins to distort the polarization 
spectra, typically ~ 150 /iT. 



IV. RESULTS 

Figure|4]shows the data obtained with the layout described 
in the last section (thick line) and the theoretical spectra (thin 
line). All experimental spectra were taken with pump and 
probe beam intensities < 0.1 7 sa t to reduce saturation effects. 
The spectra shown are of transitions from the upper hyper- 
fine level of the ground state, i.e. the 5 2 S*i/2 (F = 2) — > 
5 2 Pa/2 (F' = 1, 2, 3) transitions in 87 Rb, 5 2 S 1/2 (F = 3) -> 
5 2 P 3 /2 (F' = 2, 3, 4) transitions in 85 Rb, and 6 2 S 1/2 (F = 
4) ->■ 6 2 P 3/2 (F' = 3, 4, 5) transitions in 133 Cs. The closed 
transitions in this group are the well-known "cooling" transi- 
tion used in laser-cooling experiments. Zero detuning is rela- 
tive to the highest-frequency transition, and all detunings are 
given in units of A/2ir. The magnitude of each feature is 
given in volts, and will depend on the gain resistance in the 
photodiode circuit (1 Mil in our experiment). For Cs, two 
sets of data are shown; one taken with the vapor cell at room 
temperature (dashed line), the other with the cell in an ice bath 
as described in the previous section (solid line). 

In these spectra, polarization signals of closed transitions 
dominate, with magnitudes up to three times that of the 
next largest feature. This is true for all three species, and 
is in marked contrast with conventional saturated absorp- 
tion/hyperfine pumping spectra IU7II . As predicted, the F — > 
F' = F transitions produce a moderate amount of anisotropy, 
and the lowest frequency open transitions give the smallest 
signal. A graph showing the distribution of times of flight 
for rubidium at room temperature through a beam with the 
experimental waist shows a most probable time of flight of 
5 (is, and a significant fraction of the atoms with times of 
flight longer than 10 /is. Therefore we do indeed see that the 
timescale needed to build up anisotropy is significantly longer 
than the spontaneous decay time, and comparable to transit 
times. The accumulation of anisotropy can be accelerated 
by using a more intense pump beam; however, this will not 
greatly increase the signal as most atoms are already pumped 
into the extreme state within the transit time. This is also what 
is observed experimentally. 

Spectra taken of transitions from the lower hyperfine level 
of the ground state show a more complicated pattern. The 
(low-frequency) closed transitions in Cs and 85 Rb show a 
characteristic strong dispersion feature, but in 87 Rb the largest 
feature arises from one of the cross-over peaks (Figure|5}. In 
85 Rb closely-spaced transitions generate a polarization spec- 
trum in which individual peaks merge, making exact matching 
of features and transitions difficult. 
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FIG. 4: Experimental (thick line) and theoretical (fine line) polariza- 
tion spectra of D2 line transitions in a) 8r Rb, b) 85 Rb, and c) 133 Cs. 
Spectra were obtained by tuning the laser frequency to drive transi- 
tions from the upper hyperfine level of the atom's ground state. All 
three spectra are dominated by the strong dispersion features associ- 
ated with the closed transition. For the Cs upper hyperfine spectra, 
solid (dashed) lines represent spectra taken with the vapor cell at 0° 
(23°) C. 
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FIG. 5: Experimental (thick line) and theoretical (fine line) polariza- 
tion spectra of D2 line transitions from the lower hyperfine level of 
the ground state in 87 Rb. 



V. DISCUSSION 

By comparing theoretical and experimental traces, we see 
immediately that our model reproduces experimental features 
of upper hyperfine transition spectra with a high degree of ac- 
curacy. Fine spectral details like the "horns" resulting from 
the closely-spaced F = 2 — > F' = 2 and X3.1 cross- 
over peak in 87 Rb arise automatically from the calculated 
anisotropies. The magnitudes of theoretical spectral peaks rel- 
ative to the large closed transition peaks also agree well with 
experimental data, especially in 85 Rb. The effect of reducing 
vapor pressure by cooling the Cs cell is readily apparent. Al- 
though the central peaks are still offset from theoretical spec- 
tra compared to their counterparts in the Rb spectra, the mag- 
nitude of the offset is significantly less at °C than at room 
temperature, and the shapes of the spectra broadly agree. 

For lower hyperfine spectra (e.g. the 87 Rb spectra in Figure 
|5Jl, we observe that the sign of the anisotropy for the closed- 
transition peak is reversed compared to corresponding peaks 
in upper hyperfine spectra, as predicted in our model. How- 
ever, most other spectral features differ markedly from pre- 
dictions. Crucially, though, the model spectra show similar 
discrepancies for each of the three species, indicating that the 
lack of agreement must be due to physical processes not in- 
cluded in our model. As a test, we ran our simulation again, 
this time allowing the position and magnitude of the three 
Lorentzians (representing three transitions) to float uncon- 
strained by any input parameters. We found that no values of 
detuning or anisotropy could account for the observed spectra. 

To gain insight into the anisotropies induced by lower hy- 
perfine transitions (particularly "open" transitions, i.e. those 
with a decay channel to the other ground state level), we per- 
formed experiments designed to measure the anisotropy di- 
rectly. As in Reference |6], a quarter-wave plate was inserted 
before the PBS, oriented such that it converted the circular 
polarization component of the probe which drives <j + (<r _ ) 
transitions into vertically (horizontally) polarized light. The 
output from the arms of the PBS then directly measures the 



absorption experienced by the components driving <j + and a 
transitions. 

For upper hyperfine transitions, the largest anisotropies (i.e. 
largest difference in signal between the two arms) resulted 
from the closed transitions. This agrees with theoretical pre- 
dictions and is reflected in our data. For lower hyperfine tran- 
sitions, large and strongly negative anisotropies occurred at 
the cross-over frequencies - a fact reflected in our data, but 
not in the model. This suggests that for such transitions, it is 
no longer valid to assume that the strength of crossover fea- 
tures is the average of the two associated resonances. 

The present model assumes a uniformly distributed circular 
laser beam. The experimental beam has a Gaussian profile, 
which will lead to slight variations in the absolute height of 
the spectroscopic peaks. It will not, however, account for the 
difference between upper and lower hyperfine transitions. It is 
also possible that if the optical pumping is not complete for the 
lower hyperfine states, ignoring the coherences could have a 
large effect. The most significant approximation that is likely 
to break down for lower hyperfine transitions is the assump- 
tion that only atoms with no axial velocity will contribute to 
the anisotropy. As was shown in 11711 . a full description of the 
spectrum, especially open transitions, must take into account 
a large velocity class of atoms moving along the beam, not 
just the "stationary" ones. The contribution of these nonzero 
velocity classes will be more significant for F = T — 1/2 
because the line strength factors for the closed transition are 
weakest; for F = T + 1/2, line strength factors are strongest 
for the closed transition, and the "stationary atom" approxi- 
mation is reasonable. 

Expanding the model to include a large velocity class would 
be computationally intense, and is beyond the scope of this 
work. Walewski et al. have performed numerical simula- 
tions which included nonzero axial velocities, and success- 
fully used their model to explain features of polarization spec- 
tra in flames 1 18]. However, the organic molecules they stud- 
ied have very different structures and properties from those of 
the alkali atoms examined here, and adapting their model to 
account for alkali spectra would be correspondingly nontriv- 
ial. 

In comparing theoretical predictions with experimental sig- 
nals, we have chiefly focused on the shape and relative magni- 
tudes of spectral features. The absolute magnitude of the ex- 
perimental polarization spectroscopy signal depends not only 
on the gain resistor in the photodiode circuit, as discussed in 
the previous section, but also on the intensity of the probe 
beam in the absence of a vapor cell (see Eq. 7 in reference 
(al). This in turn depends on the optical thickness of the va- 
por, which is both temperature-dependent (as noted earlier in 
our discussion of cesium) and different for each species. Be- 
cause of these factors, the species with the largest calculated 
anisotropy will not necessarily produce the largest experimen- 
tal signal. Our decision to display theoretical and experimen- 
tal spectra on separate scales reflects these considerations. 

Finally, we note that although the peaks associated with the 
closed transitions in the upper hyperfine spectra are ideal for 
laser locking, with steep slopes centered on line center, this is 
not true for the lower hyperfine spectra, where the position of 
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each zero crossing is not trivially related to the position of the 
resonances. 



VI. CONCLUSIONS 

We have presented a model of polarization spectroscopy 
based on numerical integration of population rate equations, 
and shown how this model can be used to simulate the opti- 
cal pumping process that leads to an anisotropic population 
distribution, and thus an optically anisotropic medium. The- 
oretical polarization spectra generated by this model account 
very well for F = I + 1/2 -> F 1 transitions in 87 Rb, 85 Rb, 
and 133 Cs, but not as well for F = T - 1/2 -> F' transi- 



tions. For F =1+1/2 -> F', the closed F' = F + 1 
transition dominates the spectra, with a steep slope through 
line center which makes an ideal frequency reference for laser 
locking. Although we have not studied them experimentally 
in this work, we would expect similar results for the D2 spec- 
tra of other alkali species, e.g. 6 Li and 23 Na. 
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